Efficient DFT+f/ calculations of ballistic electron transport: Application to Au 

monatomic chains with a CO impurity 
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An efficient metliod for computing the Landauer-Biittiker conductance of an open quantum system 
within DFT+f/ is presented. The Hubbard potential is included in electronic structure and transport 
calculations as a simple renormalization of the non-local pseudopotential coefficients by restricting 
the integration for the on-site occupations within the cutoff spheres of the pseudopotential. We apply 
the methodology to the case of an Au monatomic chain in presence of a CO molecule adsorbed on it. 
We show that the Hubbard U correction removes the spurious magnetization in the pristine Au chain 
at the equilibrium spacing, as well as the unphysical contribution of d electrons to the conductance, 
resulting in a single (spin-degenerate) transmission channel and a more realistic conductance of 1 Go. 
We find that the conductance reduction due to CO adsorption is much larger for the atop site than 
for the bridge site, so that the general picture of electron transport in stretched Au chains given by 
the local density approximation remains valid at the equilibrium Au-Au spacing within DFT+U. 

PACS numbers: 73.63.-b, 73. 23. Ad, 71.28.-|-d, 72.10-d 



I. INTRODUCTION 



Electron transport through atomic-sized metallic con- 
tacts in the low-bias regime has a ballistic nature and 
is commonly studied within the Landauer-Biittiker the- 
ory. The current is carried by electronic quantum chan- 
nels, which can be partly transmitted and partly re- 
flected. The ballistic conductance of the tip-nanocontact- 
tip system is proportional to the total transmission at 
the Fermi level, G = e^/h T{Ef), where e is the elec- 
tron charge, h is Planck's constant, and T is the sum of 
the transmissions for the majority and the minority spin 
componentsiii^ For an ideal one-dimensional conductor, 
such as a monoatomic metallic chain, the ballistic con- 
ductance is proportional to the number of bands cross- 
ing the Fermi level.— When the system contains some 
source of scattering, an adsorbed impurity or structural 
disorder, for instance, the conductance is lower than 
the number of bands, so that the transmission has to 
be explicitly computed by solving an electron scattering 
problcmji In realistic systems, a convenient method for 
transmission calculations with density functional theory 
(DFT) and plane wave basis sets has been put forward 
by Choi and Ihm— in the case of Kleinman-Bylander 
pseudopotentials^ (PPs), and has been extended to ul- 
trasoft pseudopotentials^ (US-PPs) by Smogunov and 
coworkersi^ This methodology has already been applied 
to several materials, such as monatomic chains of mag- 
netic Sd-transition metals^ii tip-suspended chains of Ni^ 
Pdr- Ptr^, and Au^^ as well as to some impurity sys- 
tems such as CO on Pt monatomic chains CO on Au 
chainsH and nanocontacts^i or atomic Ni on Au chainsjii 
The more widely used nonequilibrium Green's function 
technique (see Refs. [T^ and [TgI . for instance) has been 
applied to some of these systems (e.g., tip-suspended 
chains of Aui^ and Pt^^ atop-adsorbed CO impurity on 
monatomic chains of Auii and Pi*^), usually giving re- 



sults in good agreement. 

Conventional DFT-based methods correctly describe 
the ballistic transport properties of many systems,— but 
are known to give wrong conductance values in weakly 
coupled molecular junctions (see, e.g., Refs. [l"9l - [2ll and 
references therein). This shortcoming of the DFT ap- 
proach has been recently attributed to the self-interaction 
(SI) error)^ which affects electron self-energies ob- 
tained through standard local or semi-local density 
functionalsi^ However, even in the class of systems where 
the standard exchange-correlation functionals usually 
yield good conductance values (e.g., metal nanocontacts 
and nanowires), SI errors might give rise to a wrong po- 
sitioning of the conductor electronic bands with respect 
to the Fermi level and hence alter the number of avail- 
able conductance channelsi^ In Au monatomic chains, 
for instance, the 5d-electron binding energies are lowered 
because of the SI so that d-bands are pushed toward the 
Fermi level (Ep) and two band-edges touch for Au- 
Au spacings close to the equilibrium value The spin- 
degeneracy of these bands is lifted by the Stoner instabil- 
ity induced by the extra density of states at Ep, result- 
ing in a slightly magnetic ground state and two spurious 
conductance channels in addition to the two s-channels 
(taking into account spin)i^i2i This theoretical predic- 
tion of the conductance thus gives 4e^//i for the pristine 
Au chain (or 6 e^//i in the spin-unpolarized casei^), while 
one would expect a value around Gq = 2 /h from a sin- 
gle spin-degenerate channel, as shown by experiments on 
clean Au nanocontactsi^ 

A rather simple, but very efflcient and popular way to 
tackle the SI problem is the DFT-I-C/ method, an exten- 
sion to standard DFT originally aimed at improving the 
description of electron-electron correlations of strongly- 
localized electronic states (e.g., id states in transition 
metals or 4/ states in rare earths) In this method, 
the Kohn-Sham (KS) Hamiltonian is augmented with a 
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Hubbard-C/ potential, which can be derived from a mean- 
field treatment of a many-body Hartrec-Fock Hamilto- 
nian acting on the manifold of localized orbitals. Since 
SI is absent in the Hartree-Fock method, the SI error 
ascribed to the approximate functional will be partly 
relieved for the localized electron manifold.— A previ- 
ous study on monatomic chains of 3d and Ad transition 
metals has shown that the DFT+U method can improve 
the local density approximation (LDA) description of the 
electronic structure with an accuracy comparable to that 
of a more sophisticated SI correction scheme<^ A great 
advantage oiDFT+U is that it does not add any substan- 
tial complication to standard DFT techniques and retains 
their computational efficiency in treating very large sys- 
tems, at variance with other methods more specifically 
designed to cope with SI errors, such as SI correction 
schemea^^ or hybrid functionals.— 

Some drawbacks or deficiencies, which have been only 
partly solved or addressed, can also be identified in this 
method. For instance, the choice of the parameters in 
the Hubbard Hamiltonian (just U in the simplest version, 
but many more if one considers magnetic exchange j^i^ 
spin-orbit coupling)i2i2i or inter-site interactions^) may 
critically affect the results. These parameters can be de- 
rived from rcnormalized atomic values, from constrained- 
occupation calculations (within linear muffin-tin orbital 
or similar methods2^i21i2i), or through linear- response?^ 
but in several cases one or more set of values in a rea- 
sonable range are investigated. The formulation of the 
Hamiltonian itself is subject to discussion because of 
the double counting term, which subtracts those energy 
terms that are already accounted for by the underlying 
density functionali^ The calculation of the on-site oc- 
cupations also introduces some arbitrariness: they can 
be obtained by integrating the charge inside an atomic 
sphere,— or derived from the overlap of the KS solutions 
with a set of localized wave functionsS^ (e.g., atomic-like 
states centered on the Hubbard atoms). Nevertheless, 
the method has proven useful in many cases where stan- 
dard density functionals fails and it is still widely used 
and still subject to development 

In this work, we introduce a simple and efficient 
method to calculate the Landauer-Biittiker ballistic con- 
ductance within DFT+U in a scheme based on plane 
wave basis sets and ultrasoft pseudopotentialsi^ii The 
F)FT-\-U on-site occupations arc computed using atomic- 
derived wave functions truncated at the core radius of the 
corresponding PP^^ so that the Hubbard potential can 
be rewritten in terms of the PP projectors only and it can 
be incorporated in the non-local part of the PP as a sim- 
ple renormalization of its coefficients. We shall illustrate 
the method for the widely-used rotationally-invariant for- 
mulation of the Hubbard Hamiltonianj^i^^ but a general- 
ization to more complex Hamiltonians would be straight- 
forward. 

We apply this scheme to study the effect of CO ad- 
sorption on the ballistic transport properties of an Au 
monatomic chain at the equilibrium Au-Au spacing. In- 



deed, this system was previously investigated for dif- 
ferent Au strains, but the low-strain limit could not 
be addressed because of the spurious contribution of d- 
electrons to the conductance^^ We show here that the 
Hubbard potential relieves the SI error of Sd-states in the 
Au chain, so that the corresponding bands shift toward 
higher binding energies and do not present any band-edge 
at Ep. As a consequence, a more realistic conductance 
of 1 Go and a non-magnetic ground state are recovered 
in the Au chain, thus allowing us to assess the effects of 
CO adsorption also in the low strain limit. We find that 
the conductance reduction is substantially larger for the 
atop adsorption than for the bridge adsorption because 
of the different position of the transmission dip due to the 
hybridization resonance between the Scr-molecular level 
and the s-band of Au, much closer to Ep in the atop ge- 
ometry. This confirms the general picture given by the 
LDA study at larger strains 

The paper is organized as follows: In Sec. [Ill "we provide 
the theoretical basis for embedding the Hubbard poten- 
tial in the electronic structure and transport calculations 
through the pseudopotential coefficients. In Sec. IIIIl we 
apply the method to the electronic, magnetic, and struc- 
tural properties of the pristine Au chain. Then, in Sec. lIVI 
we report the LDA+U adsorption energetics of CO on the 
Au chain and we investigate the effect of the Hubbard U 
on the transport properties. Finally, our conclusions will 
follow in Sec. El 

II. METHOD 

A. A scheme for ballistic transport within DFT4-I7 

Since its earlier versions j^i^i^ several variants of the 
DFT+C/ method have been proposedi^i^i^ In order 
to illustrate our method, here we cho ose the simplifi ed 
rotationally-invariant formulation of jPudarev et a/.l .— 
which is widely used and has already been imple- 
mented in a plane wave-pseudopotential codci^i^ In the 
rotationally-invariant Hubbard Hamiltonian, the main 
effect of the on-site Coulomb repulsion is described 
through a single parameter U and the Hubbard energy 
term that is added to the usual DFT total energy reads 
as 

/,cr mm' 

where n^^^^, are the on-site occupation matrices for 
atomic site / and states of spin cr. The indexes m and m' 
run from —I to where the orbital angular momentum 
I is fixed by the choice of the localized electron manifold 
{1 = 2 for the d shell, / = 3 for the / shell, . . . ) and is 
omitted here for shortness of notation. The occupation 
matrices can be written in a quite general form as: 

^mm' = E (^L I Prnm' I K.) , (2) 



3 



where ip^^ are the solutions of the KS equation with 
spin (7, /k« are single-particle occupation factors, and 
P,nm' generahzed projection operators.— In plane 
wave codes, the latter can be conveniently chosen as fully 
separable, 

PLn'^\vL){V>Ll (3) 

where the (pl^ can be, for instance, atomic pseudo-wave 
functions. Consequently, the Hubbard potential 

/ J / J * mm' mm' i 
I mm' 

will appear in the KS equation, with V^l^^, defined as 

^mm' ^ ^ I,<j ^ 2 v"™'™ ^'^m'mJ- \P ) 

'^'"'mm' 

From Eq. Q we notice that the Hubbard potential has 
the same form as the non-local part of a pscudopotential, 
which in the ultrasoft case is usually written as^ 

/ I] 

where -D^j'^ are the screened coefhcients and /?/ are the 
projector functions associated to atom /. However, while 



PP projectors are nonzero only within a sphere of radius 
around the atom, the DFT+U projectors in Eq. ((S]) 
have in principle an infinite support, even if in practice 
the atomic wavefunctions (r| (^^J decay rapidly to zero 
sufficiently far away from the atomic center /. This fea- 
ture of the Hubbard potential is not desirable because it 
makes the solution of the scattering problem more com- 
plicated and it would also increase its computational size. 



This kind of problem can be overcome by limiting 
the integration for the on-site occupation matrices 
within a sphere of finite radius centered on atom /, thus 
neglecting the contributions coming from the charge out- 
side the spheres. A similar approach has been proposed 
in the projecto r augmented-wave (PAW) framework by 
iBengone et al\ ^ who showed that the excluded charge 
is usually small for the localized states in the Hubbard 
Hamiltonian and that it is possible to use the PAW pro- 
jectors instead of atomic wave functions in the Hubbard 
potential by adopting this approximation. The similar- 
ity between the US-PP method and the PAW formal- 
ism allows us to apply a PAW-like transformation^ to 
the all-electron (AE) version of our projection operators, 

^mm' — \ 'Pm'^){4'm'^ I: Order to obtain the correspond- 
ing pseudo (PS) version^^ as in Ref. HI: 



m Pmm' I ^k^) = ('^k.l P:km' I V'k.) + E^V'k.l A) [('/'f I P:km' I - ('An P^"m' I C)] i^AM (7) 

ij 



where (f>^^ i(f>T) the AE (PS) partial wave correspond- 
ing to the PP projector /3i of atom /, and V'ku are the PS 
wavefunctions that are obtained from the solution of the 
KS equation with US-PPs (we have omitted the atom in- 
dex / for shortness of notation) . The integrals needed for 
the AE and PS matrix elements between square brack- 
ets can be performed within the augmentation spheres 
around the selected atom (|r — R/| < Tc), since AE and 
PS partial waves coincide outside the spheres by con- 
struction. 

If we assume that the projection operators P^^, are 
sufficiently localized within the atomic spheres and that 
partial waves and projectors form a complete basis in- 
side those regions, we can apply to P,^"^; the following 
equality: 



= (^k. I P I ^k«) - ^ (V'k. I ft ) (C I B I 0f ) (/3j I ^k„) , 

ij 

(8) 

which strictly holds for any arbitrary operator B entirely 
localized within the atomic spheres. This allows us to 
obtain from Eq. 1^ an approximate expression for the 



projection operators: 

iKllPmm' I Kl) ^ E(V^k.| ft)(0f I PZ.' I 0r)(/3.l ^k.) 
ij 

(9) 

where we neglected the contribution of the atomic wave 
functions outside the augmentation spheres. If the above 
expression is used in Eq. (|4]), the Hubbard potential can 
be rewritten in the following form: 



EE 

I ij 



^mm' Wi I Pmm' I </'i ) 



(10) 

and it can be readily incorporated in the non-local part 
of the US-PP, resulting in 



NL 



I ij 



(11) 



where A[j'^ are the quantities between square brackets 
in Eq. (fTn|) . The effect of the Hubbard potential can 
thus be included in the KS equation through a simple 
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rcnormalization of the non-local PP cocfScicnts and no 
additional projectors other than those already required 
by the US-PP are needed. The equation for the electron- 
scattering problem in the framework of US-PP and plane 
waves can be written as (Rydberg atomic units, e^/2 = 
2m = h = 1, are used)!^'^'^^ 



[- V2 + T/cff + V^NL]l*k) =S|*k) , 



(12) 



transmission is calculated only for the "pseudo" method 
and the ballistic conductance has been obtained within 
the Landauer-Biittiker formalism by evaluating the total 
transmission at the Fermi energy, G = e^/h ■ T{Ep). We 
do not adopt here any self-consistent determination of 
U but we rather compute the electronic structure for 
some values of t/ in a range of interest. 



where obtained from in Eq. ^ by replac- 



ing the screened coefficients with D 



D 



(the qfj 's being the integrals of the augmentation func- 
tions defined in Ref. Q). Therefore, the only additional 
step to include the Hubbard potential in this ballistic 
transport scheme is to replace in Eq. (fT2|) with a 
new potential V^l+u where the coefficients are given by 
Dlf + A^f -Eqlj. In this way, the T)¥T+U transport 
calculation do not present any additional theoretical or 
technical difficulty than those already discussed and re- 
solved in Refs. iand0. 



B. Computational details 

We have implemented the DFT-f-C/ transport method 
outlined above in the Quantum ESPRESSO package,— 
integrating it with the existing implementations of the 
DFT+J7 methodSS, (in the PWscf code) and of the baUis- 
tic transport with US-PP (PWcond code^). The DFT-t-fJ 
transmission calculations proceed in two steps: first, the 
ground-state electronic structure of the system is com- 
puted with PWscf to obtain the local potential VcS, the 
screened coefficients dI'^ ^ as well as the Hubbard coef- 
ficients A,[j'^. In this step, the A,[j'^ must be updated at 
each iteration of the self-consistent loop, similarly to the 
US-PP screened coefficients, because of their dependency 
on the on-site occupations. Second, the potential Vp^^^u 
is used to solve the scattering problem [Eq. ((T2|) ] with 
the same techniques of Rcf. 7. 

In Sec. IIIH we will apply this method to the isolated 
Au monatomic chain, while in Sec. IIVI we will study the 
ballistic transport through the chain in presence of a CO 
impurity adsorbed on it. The calculations have been per- 
formed within the LDA using the same computational 
parameters described in Ref . [l3l for the CO/Au chain sys- 
tem. In particular, we use here the same US-PP of Au, C, 
and O, as well as the same plane wave cutoffs for the wave 
functions and the charge density presented there. For the 
spin-polarized calculations of the infinite Au chain we re- 
duced the smearing parameter for the electronic occupa- 
tions to 0.005 Ry and sampled the irreducible Brillouin 
zone with 46 fc-points. In our LDA-I-J7 calculations, the 
Hubbard potential is applied to the 5d-electron manifold 
of Au, either using the full atomic pseudo- wave functions 
as in Eq. ([5]), or using the PP projectors as in Eq. (fTT|) . 
For convenience, we named here "atomic" method the 
former and "pseudo" method the latter. The electron 



III. ELECTRONIC AND MAGNETIC 
PROPERTIES OF THE AU CHAIN 

In this section we discuss the effects of the Hubbard po- 
tential on the structural, magnetic, and electronic prop- 
erties of the infinite Au monatomic chain computed in 
the LDA. We will compare the plain LDA results with 
the LDA+J7 results obtained with both the atomic and 
pseudo methods (see above) for the on-site occupations. 
Finally, we will validate the complex band structure 
(CBS) calculation for the pseudo LDA+[/ method by 
showing that the usual band structure obtained with pe- 
riodic boundary conditions (PBC) along z is reproduced 
by the CBS at real fc^-points. 
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FIG. 1. (Color online) Cohesive energy (a) and magnetic mo- 
ment per atom (b) for the pristine Au chain as a function of 
the Au-Au spacing, dchain. The LDA-I-C/ results obtained with 
three different values of f/ (1.0 eV, 3.0 eV, and 5.0 eV) are pre- 
sented for each of the two projection methods ( "atomic" and 
"pseudo", see text). The corresponding equilibrium Au-Au 
spacing dcq (in A) is reported in the insets of panel (a). In 
(b), the magnetic moments obtained from plain LDA calcu- 
lations are also reported for comparison (black diamonds). 
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FIG. 2. (Color online) Electronic band structure of the Au 
monatomic chain at the LDA equilibrium spacing for differ- 
ent values of the U parameter (in eV) and projection method 
("atomic" or "pseudo"). In the leftmost panel the band struc- 
ture from plain LDA is reported for comparison, while in the 
rightmost panel different values of U have been used for the 
two projection methods. The bands are plotted for fcz in- 
side the irreducible Brillouin zone (0 kz < n/a^, with 

da = dchain). 



In Fig. [TJ we report the ground-state energy and mag- 
netic moment per atom for the pristine infinite chain 
as a function of the Au-Au spacing, dchain- We study 
the dependence of these quantities on the strength of 
the Hubbard potential for several U values up to 5 eV. 
For each value of U, the total energy of the isolated Au 
atom computed consistently is used as reference energy, 
hence the plots in Fig. [T]Ja) represent the cohesive en- 
ergy of the chain. As U is increased, the equilibrium 
spacing of the Au chain (see deq in the insets) slightly 
expands from the plain LDA value (2.51 A) until it 
saturates to about 2.55 A (2.54 A) in the atomic (pseudo) 
method for U ^ 3.0 eV. The cohesive energy instead de- 
creases for both projection methods, but more rapidly 
for the pseudo method than for the atomic one. The 
total ground state magnetic moment per atom for spac- 
ings 1.9 A ^ rfchain ^ 2.8 A is reported in Fig. [T]Jb) for 
both pseudo and atomic methods. The plain LDA result 
(diamonds) shows that for 2.22 A < dchain < 2.78 A the 
ground state of the Au chain bears a finite magnetiza- 
tion. The magnetic moment per atom reaches a max- 
imum value of 0.23 /xb at dchain = 2.3 A, while at the 
equilibrium spacing it is about 30% smaller (0.15 /.te). 
When the Hubbard interaction is turned on, we observe 
a significant shift of the magnetic instability region to- 
ward smaller Au-Au spacings and a gradual reduction of 
the maximum magnetic moment as U is increased. Small 
quantitative differences between the two methods can be 
recognized, such as a more pronounced contraction of the 
magnetic instability region and of the maximum mag- 
netic moment in the atomic method with respect to the 
pseudo method. However, for strong enough Hubbard 
potentials (U ^ 3eV), both atomic and pseudo methods 
completely suppress the spurious magnetization of the 
Au chain at its equilibrium spacing. 
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FIG. 3. (Color online) Complex band structure of the Au 
monatomic chain within LDA-t-t/ {U — 3eV) from PWcond 
(circles and points) and band structure with full PBC from 
PWscf (solid lines). The central panel contains the CBS eigen- 
values at real fcz (circles), while those at imaginary and at 
kz = 1/2 + iK are on the left and right panels (points), re- 
spectively (see also Ref. 0) ■ The wave vector kz is expressed 
in units of 27r/a3, with as — dchain ~ 2.51 A. 



In Fig. [21 we report the spin unpolarized electronic 
band structure of the Au chain obtained with plain LDA 
and with LDA-I-J7 for selected values of U. The large 
electron density of states at Ep, which is the main rea- 
son for the magnetic instability seen in the LDA, is due 
to the upper edge of the doubly-degenerate d^z/dyz band 
pinned at Ep. Within LDA+U, these bands, as well as 
the other filled bands (mainly of d-character) , are pro- 
gressively pushed toward higher binding energies as U is 
increased, resulting in the suppression of the magnetic 
instability and of the spurious d-channels for U values 
as small as 3eV. This effect is qualitatively the same 
for both atomic and pseudo methods, but the downward 
shift of the d-bands is more pronounced when using the 
full atomic wave functions4S A perfect matching of the 
"atomic" and "pseudo" band structures is not possible 
because of the different broadening of the bandwidths, 
but a comparable energy shift can be obtained if a larger 
value of U is used in the pseudo method (see rightmost 
panel of Fig. [2]). 

In order to perform ballistic transport calculations 
within DFT+U, the CBS of the leads has to be com- 
puted at the same level of accuracy as the usual elec- 
tronic band structure calculated using full PBC. We test 
if our implementation meets this requirement by comput- 
ing the CBS of the Au monatomic chain within LDA-I-J7 
for a representative value of [/ = 3 eV. In Fig.[3l the CBS 
computed with the ballistic transport code PWcond (sym- 
bols) is compared with the band structure obtained with 
the electronic structure code PWscf (full lines). On the 
energy scale used in the figure, the CBS at real (cen- 
tral panel) is indistinguishable from the band structure 
with full PBC. Therefore, when used in transport calcu- 
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TABLE I. Optimized distances (in A) and chemisorption energies (in eV) for the bridge and atop adsorption geometries 
within LDA and LDA+C/. The labels "atomic" and "pseudo" refer to on-site occupations computed using, respectively, atomic 
wavefunctions and PP projectors. The reference energy of the isolated chain is consistently computed within LDA+(7, using 
the same value of U and type of projectors. 



Bridge Atop 
atomic pseudo atomic pseudo 



u 


rfAu-C 


dc-o 


-£'chcm 


dAu-C 


rfc-o 




dAu-C 


dc-o 




dAu-C 


dc-o 




0.0^ 


1.952 


1.156 


-2.4 








1.891 


1.136 


-1.0 








1.0 


1.957 


1.155 


-2.4 


1.952 


1.156 


-2.4 


1.901 


1.135 


-1.0 


1.896 


1.135 


-1.0 


3.0 


1.971 


1.151 


-2.1 


1.953 


1.154 


-2.1 


1.931 


1.133 


-0.8 


1.920 


1.134 


-0.8 


5.0 


1.988 


1.148 


-1.7 


1.959 


1.151 


-1.6 
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LDA results from ReL [H. 



lations this CBS will result in a single (spin-degenerate) 
transmission channel open at for the Au chain, in 
more reasonable agreement with what is observed in ex- 
periments. 



IV. AU CHAIN WITH CO IMPURITY 
A. Geometry and energetics 

We started by computing the optimized structures of 
the bridge and atop adsorption geometries of CO within 
LDA+C/ for several values of U , following the same cri- 
teria described in Ref. [l^ The molecule is kept in an 
upright position and the positions of C and O are op- 
timized, while all Au atoms are kept fixed and aligned, 
equally spaced by the LDA equilibrium distance. 

In Table m we report the optimized C-Au and C-0 dis- 
tances (dAu-c and cZc-O ; respectively) and the chemisorp- 
tion energy (E'chem) for the two geometries. These re- 
sults present a general trend of a progressive weakening 
of the interaction between the CO and the Au chain for 
increasing U values in both geometries, as shown by the 
significant decrease of i?chcm, the increase of dAu-c, and 
the shortening of o?c-o for U > 1 eV. Notice that the 
lowering of -BchGin due to the Hubbard U correction docs 
not change the energetic preference for the bridge site 
that was found in previous calculations^^ This conclu- 
sion does not depend on the projection method used, 
because the differences in £^chom between the atomic and 
pseudo methods stay below 0.1 eV, about one order of 
magnitude smaller than the bridge-atop -Bchcm difference 
(~leV). 

The weakening of the CO-Au interaction can be ratio- 
nalized in terms of the energy shift of the Au hd bands: 
larger values of U will result into greater shifts of the 
d-band center toward higher binding energies, further 
and further away from the Fermi level. This in turn 
will decrease the reactivity of the metal toward adsor- 
bates, resulting in lower adsorption energies, larger ad- 
sorbate/metal binding distances, and smaller molecular 
bond relaxationsi^ 



B. Ballistic conductance 

In the previous section, we have shown that LDA+C/ 
gives the correct number and type of conductance chan- 
nels for the Au monatomic chain. Here we will use this 
improved description of the CBS of the Au chain to study 
the effects of CO adsorption on the ballistic electron 
transport through the chain at low strains (dchain — 2.5 A 
within LDA). We shall first present the LDA+C/ trans- 
mission for a selected value of U and compare it with the 
plain LDA result then we will study the dependence 
of the tipless ballistic conductance on U . The bridge and 
atop geometries obtained previously with plain LDAji^ 
will be used here, because we want to discuss the ef- 
fect of U through the electronic structure, rather than 
through the changes of the adsorption geometry (which 
arc nevertheless small, sec Table |T|. 

In Fig. [31 the transmission through the Au chain with 
the CO impurity and the number of available channels 
(i.e., the transmission of the pristine infinite chain), are 
reported as a function of the electron scattering energy. 
The figure displays both the LDA+C/ results obtained 
with U = 3eV and the previously obtained LDA re- 
sults, reported here for comparison. The removal of the 
spurious conductance channels due to the dxz/dyz bands 
(cf. Fig. [2]) gives rises to important changes in the ballistic 
conductance, which becomes 1 Go for the pristine chain 
within LDA+C/, and drops below 1 Gq in presence of CO. 
The transmission of the remaining spin-degenerate s- 
channcl is affected by CO in a rather different way in the 
two adsorption geometries and results in a much larger 
conductance for the bridge geometry (G ~ 0.7 Go) com- 
pared to the atop geometry (G ~ 0.3 Go). The smaller 
conductance in the atop geometry stems from the dip 
in the s-transmission associated to the 5iTa-hybridization 
resonance which, at variance with the bridge geometry, 
falls very close to E-p. This feature is also found in the 
LDA transmission, but the conductance-cutting effect of 
this interaction is hidden by the spurious conductance 
due to the extra d-channels at E-p. 

At other scattering energies we also notice some 
changes in the transmission, especially below E-p^ where 
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FIG. 4. (Color online) Electron transmission (solid lines) and number of available transmission channels (dashed lines) as a 
function of the scattering energy for the bridge and atop geometries (left and right panels, respectively). The figure shows both 
the LDA+(7 result with U = 3.0 eV (thick red lines) and the plain LDA resuHii^ (thin blue lines). 



the d-channcls give the largest contribution. In general, 
the LDA+C/ transmission is slightly lower than the LDA 
one, excepting some energy regions (e.g. around —3.5 eV) 
where the number of channels increases in the LDA+J7 
because of the d-band shift to higher binding energies 
induced by the Hubbard potential. Below — 3.5eV, a 
single s channel is present and the main differences are 
localized in correspondence of transmission dips or at the 
lower band edge. Above E-p, the number of transmission 
channels is not influenced by the Hubbard potential, nev- 
ertheless in the bridge geometry the transmission is mod- 
ified between 2 eV and 3 eV because of a transmission dip 
shifting by about 0.2 cV toward Ep. This dip is due to the 
antibonding 27r* resonance^^^ which is more coupled to 
the metal d-states with respect to the 5(Ta resonance and 
is thus sensitive to the energy shifts of the d-bands. 

Finally, we discuss how the tipless ballistic conduc- 
tance of the bridge and atop geometries depends on the 
strength of the Hubbard potential. In Fig. [Sj we report 
the conductance for some selected values of U between 
OeV (plain LDA) and 5eV. For small values of U, the 
tipless conductance decreases rather rapidly with U be- 
cause of the vanishing of the d-channel contributions to 
the conductance, and already at U = I both geometries 
have a conductance lower than 1 Gq. For U ^ leV, the 
conductance as a function of U stays almost constant 
in the bridge geometry, while it further decreases in the 
atop geometry as larger values of U are considered. How- 
ever, for [/ > 3 eV the two geometries have well-separated 
conductances, irrespectively of the value of U. We can 
therefore conclude that the ballistic conductance in the 
bridge geometry is substantially higher than in the atop 
geometry not only at high strains, but also when the Au 
chain is not stretched. 



V. CONCLUSIONS 

In conclusion, we have addressed the self-interaction 
problem in the electronic structure and in the ballistic 
transport through a simplified DFT-I-J7 method, which 
leads to an efficient calculation of the electron trans- 





Bridge : 
Atop — e— : 







1 2 3 4 5 

Hubbard U (eV) 

FIG. 5. (Color online) Ballistic conductance in the bridge 
(blue points) and in the atop geometries (red circles) obtained 
for a selected set of U values (lines connecting the calculated 
conductances have been drawn to guide the eye). 



mission coefficients. Using an approximate value for the 
on-site occupations entering the Hubbard Hamiltonian, 
we could include the effect of the Hubbard potential in 
the Kohn-Sham equation and in the scattering equation 
through a simple renormalization of the coefficients of 
the nonlocal pseudopotential. 

We showed that the method is capable of suppressing 
the spurious magnetic instability of an Au monatomic 
chain arising at the equilibrium spacing in the LDA, thus 
recovering a non-magnetic ground state and removing at 
the same time the spurious contribution of Au d-states 
to the conductance. We find that the chemisorption en- 
ergies of CO on the Au chain can be substantially lower 
in the LDA+U compared to the LDA, but the energetic 
preference for the bridge site previously obtained^^ with 
plain LDA is confirmed by the LDA+U calculations. The 
comparison of the ballistic conductance of the bridge and 
atop adsorption geometries, which is not feasible in the 
LDA without stretching the chain, reveals that the con- 
ductance cut mechanism seen at high Au strain takes 
place also at low strain. Indeed, in the atop geometry 
the s-transmission dip due to the 5aa resonance is very 
close to Ep and is responsible for the large conductance 
reduction, while in the bridge geometry the dip falls at 
higher energies and affects the conductance to a much 
smaller extent. 
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